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Abstract 

A one-dimensional model of a dispersive medium with intrinsic loss, compen- 
sated by a parametric drive, is proposed. It is a combination of the well- 
known parametrically driven nonlinear Schrodinger (NLS) and complex cubic 
Ginzburg-Landau equations, and has various physical applications (in partic- 
ular, to optical systems). For the case when the zero background is stable, 
we elaborate an analytical approximation for solitary-pulse (SP) states. The 
analytical results are found to be in good agreement with numerical findings. 
Unlike the driven NLS equation, in the present model SPs feature a nontrivial 
phase structure. Combining the analytical and numerical methods, we identify 
a stability region for the SP solutions in the model's parameter space. Gener- 
ally, the increase of the diffusion and nonlinear-loss parameters, which differ the 
present model from its driven-NLS counterpart, lead to shrinkage of the stabil- 
ity domain. At one border of the stability region, the SP is destabilized by the 
Hopf bifurcation, which converts it into a localized breather. Subsequent period 
doublings make internal vibrations of the breather chaotic. In the case when 
the zero background is unstable, hence SPs are irrelevant, we construct station- 
ary periodic solutions, for which a very accurate analytical approximation is 
developed too. Stability of the periodic waves is tested by direct simulations. 
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1 Introduction 



One of paradigm models adopted for the description of pulse dynamics in non- 
linear dispersive media with intrinsic loss is the parametrically driven damped 
nonlinear Schrodinger (NLS) equation. In the one-dimensional case, it is 

iUt + ^U^j:x + ^\u\'^u ^ {k - i)u + ^u* , (1) 

where u(x,t) is the complex wave field, the asterisk stands for the complex 
conjugation, real k is the frequency mismatch of the parametric drive, real 7 is 
the drive's amplitude, which may be assumed positive, and the loss parameter, 
as well as the ones accounting for the spatial diffraction/dispersion and nonlinear 
frequency shift, are normalized to be 1. Equation finds applications to 
various physical problems, such as dynamics of spin waves in magnetics 
plasma waves j2j: and Faraday surface waves in a fluid layer driven by a time- 
periodic external force Arguably, the most realistic application of eq. ^ 
is to models of nonlinear optics, such as a cavity filled with a Kerr-nonlinear 
lossy material, the driving term being provided for by an external field supplied 
at the frequency coinciding with that of the signal field, see, e.g., ref. 4 and 
references therein. 

It is well known that eq. JP) gives rise to two exact solutions for phase-locked 
solitons|S], 

u{x) = A{x)e'\ A{x) = 2y^kTTcos20sech (^a/2(A: 7cos20)a;) , (2) 

where the phase constant takes values 

Oi = sin-i (I/7) , 02 = i [tt + sin-i (I/7)] . (3) 

The solitons exist provided that 7 > 1. 

The soliton with 9 = 62, corresponding to the smaller value of the amplitude, 
^0 = 2-yAr+7^cos2^, is always unstable, while the other one may be stable. The 
soliton stability in this model was studied in detail (6j. In particular, an obvious 
necessary condition is stability of the zero solution, which implies 

k > v/72 - 1. (4) 

Thus, the stable soliton may exist in the interval 1 < 7^ < 1 -I- /c^ of the values 
of gain. Actually, the soliton is stable only in a part of this region, while in 
another part it is destabilized via a Hopf bifurcation, which was also studied in 
detail, see an earlier work and a systematic analysis in refs. 8. 

Another important dynamical model for pattern formation in nonlinear lossy 
media is based on the complex cubic Ginzburg-Landau (CCGL) equation, which 
may be written as 

iut + {Pi - il32)uxx + (ci -I- ic2)\u\^u = {k + i)u, (5) 
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where the coefficients /3i_2 and ci^2 are real, and /92 and ci cannot be negative. 
The hnear-gain coefficient on the right-hand side of eq. Q is normahzed to be 
1 [cf. the normahzed hnear-loss coefficient in eq. (Q]. The CCGL equation 
has many apphcations, including description of the onset of turbulence in the 
Poiseuille flow and Langmuir waves in plasmas ^U]. Additionally, if t and 
X are replaced, respectively, by the propagation distance z and reduced time 
r, eq. lO describes the propagation of light signals in a pumped nonlinear 
optical waveguide, with the dispersion coefficient /3i, Kerr coefficient ci, filtering 
coefficient ^i , and coefficient of two-photon absorption C2 ^ . 

Various solutions to eq. © and their physical applications have been studied 
in detail, see a recent review 11^ . A particular solution, which is unstable but 
nevertheless quite meaningful (see, e.g., a detailed consideration in ref. 13), is a 
solitary pulse (SP) that can be found in an exact analytical form |^ ^] 

u — Aq [sech (kx)]^^*^ exp (—iutt) , (6) 

where all the real constants Ao,K,fi and lo are uniquely expressed in terms of 
parameters of eq. JSJ. It is noteworthy that the NLS equation proper, and its 
damped parametrically driven version admit soliton solutions only under 
the condition which, in terms of eq. lO, is /3iCi > (in terms of nonlinear optics, 
it implies that the nonlinear waveguide is either self-focusing with anomalous 
group- velocity dispersion, or self-defocusing with normal dispersion). However, 
except for the case /32 = C2 = 0, the general CCGL equation (O has the exact 
SP solution in the form ^ irrespective of the sign of the product (3iCi, the only 
additional condition being that (3i and ci must not simultaneously vanish. The 
most important difference of the SP solution ^ from the phase-locked soliton 
(|2Jl is the presence of the chirp /i, which accounts for a nontrivial intrinsic phase 
structure in the SP. 

A natural generalization of both models and (0) is a parametrically driven 
CCGL equation, which we take in the form 

iut + iPi - i(32)uxx + (ci + ic2)\u\'^u = {k ~ i)u -f ^u* . (7) 

It may apply to the above-mentioned physical systems under properly modi- 
fied conditions. For instance, in the case of the nonlinear optical waveguide, 
eq. lO describes a situation when the loss is compensated not by intrinsic 
gain, but rather by a co-propagating pump wave. It should be mentioned 
that a parametrically driven model of that type was already introduced in a 
two-dimensional situation, which corresponds to the light propagation in a pla- 
nar nonlinear waveguide |14| . and stable anisotropic two-dimensional SPs were 
found. In one dimension, a CCGL system with a nonlocal parametric drive was 
considered recently jl5| . but the straightforward one-dimensional model (0 has 
not been considered yet. Complex dynamics of interfaces was also studied using 
this type equation |16) . 

The existence of the stable soliton ^ in the driven-NLS limit of the model 
((TJ, and of the unstable SP © in its CCGL counterpart, suggest a possibility 
to find a SP solution to eq. Q and investigate its stability, which is the subject 
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of the present work. One may expect that an essential difference of the SP 
solution from its NLS phase-locked counterpart ((^J is a nontrivial intrinsic phase 
structure (chirp). 

In section 2, we report results obtained for SPs in the model Q. First, we 
develop an analytical approximation, that predicts all the properties of SPs, 
except for the chirp, with good accuracy. Then, we demonstrate numerical 
results, which are summarized in the form of stability diagrams in the model's 
phase space. The stability domain shrinks with the increase of the diffusion and 
nonlinear-loss parameters /3i and ci in eq. {Tj) , which make it different from the 
NLS counterpart JQ). At one of the domain's borders, the soliton is destabilized 
via a Hopf bifurcation, which gives rise to a stable breather. The breathers are 
also briefly considered in section 2. 

In the case when the condition for the stability of the zero solution is 
not met, SPs keep to exist, but they cannot be stable. In this case, relevant 
solutions are periodic waves. We briefly consider them in section 3, also first in 
an approximate analytical form, and then by means of numerical methods. In 
this case, analytical and numerical results are virtually indistinguishable. 

2 Solitary-pulses 

2.1 Analytical approximation 

To develop an analytical approximation to SP solutions, we fix two parameters 
in eq. (0), /3i — ci = 1/2, which can always be done by means of rescaling 
(provided that j5\ and c\ are positive). We then assume an approximate form 
for the solution which mimics the exact solution l^Jl available in the NLS limit 

o, 

u{x) = A{x)e'\ (8) 

with constant 6. Then, substituting the ansatz (jSJ into eq. we formally 
obtain two different equations for the single real function A{x): 

^ =2[/c + 7COs(20)]^-A3, (9) 

/32^ = [l+7sin(20)]^ + C2A3. (10) 

We consider the case that /32, 1 + 7sin(20) and C2 are small. Eliminating the 
derivative term from eqs.© and (|10|l . we arrive at a simpler relation, 

(C2 + 132) A"" = (2/32 [fc + 7 cos {26)] - [1 + 7 sin(20)]) A. (11) 

An exact solution to eq. lO is given by the expression for A{x) from eq. (0), 
which we will use here, but without adopting the relation (PJ. Instead, we 
substitute the expression for A{x) into eq. Hll|l. then multiply all the equation 
again by A{x), and integrate the result over x from —00 to -l-oo. This procedure 
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leads to the following relation, that we use as a definition of 9 in the waveform 
which is now considered as an approximate solution to eq. (|7|): 

2(fc + 7COs26l)(4c2+/32) + 3[l + 7sin(26l)] = 0. (12) 

In particular, eq. H12|l . with regard to eq. Q, predicts that a stable SP solution 
may exist in the interval 



V7 l<A:<fc„_ 2(4c2+/52) ' ^ ^ 

With the increase of the combination (4c2 +132), fccr attains a minimum value. 



= v/t^, (14) 

at 

4c2 + /?2 = (3/2) Vt'-I (15) 

[at this point, the interval (|13|l collapses into a point]. With the further increase 
of (4c2 + /32), kcr grows up to a maximum value, (fccr)max = ^cr (4c2 + /32 = cxd) = 
7- 




Figure 1: A typical example of the stable solitary pulse. The panels (a) and 
(b) show the amplitude and phase profiles of the stationary solution with 7 = 
1.5, fc = 2 and P2 = C2 = 0.02. In (a), the solid and dashed curves represent 
the numerical result and analytical approximation [based on Eqs. lO and H12|l]. 
respectively. 



2.2 Numerical results 

The stationary version of eq. 10 was solved numerically, and the findings were 
compared to those obtained in the approximate analytical form based on eqs. (|2J) 
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Figure 2: The amplitude of the solitary pulse (a) and its average phase (b), as 
defined in Eq. (|15|l . vs. the parameter P2 = C2. Other coefiicients in Eq. (7J) 
are fixed: 7 = 1.5, k = 2. The symbols show numerical results, while continuous 
curves are predicted by the analytical approximation. 



and H12|l . A typical example is displayed in Fig. 1, where the panel (a) shows the 
numerically found profile of \u{x)\ together with the analytical approximation 
to it, and the panel (b) is the corresponding phase profile, which was found 
only in the numerical form, as the analytical approximation assumes a constant 
phase. 

Results for the stationary SP solutions can be summarized in the form of the 
plot displayed in Fig. 2. Its panel (a) shows the soliton's amplitude, as found 
numerically and in the analytical approximation, vs. the coefficient P2 = C2, 
which measures the deviation of the model O from the NLS counterpart (Q, 
and the panel (b) shows the average value of the SP's phase (pix), which is 
defined as 

^^ c::mHx)\'dx ^ ^^^^ 

I-^ \u{x)\^dx 

in comparison with the analytical prediction for the average phase, which is 
realized as the phase 9 defined by eq. H12|l . The numerical and analytical results 
displayed in Fig. 2 terminate at points where, respectively, the numerical solu- 
tion and eq. H12|) cease to produce any nontrivial result (the numerical solution 
shows that, beyond the termination point, only the trivial state, u = 0, can be 
found) . 

The results are further summarized in the diagram shown in Fig. 3, which 
displays a region in the parameter plane (/32 = C2, A:), in which stable stationary 
SP solutions were found [the stability was verified by direct simulations of the 
full time-dependent equation The upper right border marked by stars 
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Figure 3: The region in the parametric plane with the coordinates [32 = C2 and 
fc (7 = 1.5 is fixed) in which stable solitary-pulse solutions were found. The 
curves connecting symbols which mark stability/existence borders are guides 
for the eye. The dashed curve, close to the upper right border, is the analytical 
prediction H13() . 



is the one beyond which only the zero solution can be found numerically; the 
dashed curve is an analytical prediction for the same border as given by eq. (|13() , 
i.e., the boundary beyond which eq. H12(l yields no solutions. Below the lower 
border, which is marked by crosses, stable SPs cannot be found. The width 
of an initial pulse increases and the pulse spreads out to a spatially extended 
state. (The spatially extended state is not a uniform state.) The lower border 
is slightly upper than the instability line k — \/ ^(^ — \ of the zero solution. 
The stability region in Fig. 3 ends approximately at (/32)cutoff' ~ 0.16, where the 
upper border approaches the lower border. Lastly, the upper left border, marked 
by rhombuses, is a line of the Hopf bifurcation destabilizing the SP: behind this 
border, a stable solution is a regularly or chaotically vibrating breather, see 
below. 

The consideration of the above results, and the comparison with the ana- 
lytical prediction shows that, first, the analytical approximation turns out to 
be fairly accurate, and, second, the stability region quickly shrinks with the 
increase of the diffusion and nonlinear-loss coefficients (ii and ci in eq. lO. The 
stability region in Fig. 3 greatly extends in the direction of larger values of k at 
small finite values of (ii — ci [which may be easily explained by the analytical 
result showing that the maximum value of fc, as given by eq. I|13|) . diverges as 

/32 = C2 ^ 0]. 

Another relevant representation of the SP stability region is in terms of the 
parametric plane fc) for fixed values ii2 — C2 — 0.025, c\ — 0.5 and 7 — 1.5 of 
the other parameters, which is shown in Fig. 4. The two borders in Fig. 4 denote 
the same types of bifurcation curves as in Fig. 3. The upper border marked by 
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Figure 4: The solitary-pulse stability region in the parameter plane (/3i,fc) at 
fixed values ci = 0.5, /32 = C2 = 0.025, and 7 = 1.5. The pulses are stable 
between the two borders. The dashed curve, close to the upper right border, is 
the analytical prediction. 



stars is the one beyond which only the zero solution can be found numerically; 
the dashed curve is an analytical prediction fccr = + (8c2 + ^2/ (iiY ~ 

3}/(8c2 + /32//?i), which is a modified one of eq. The analytical prediction 

is good for large however, the approximation breaks down for small (3i, 
since the derivation of Eq. is based on the conditions (32 Pi and C2 ^ Ci . 
Below the lower border, which is marked by crosses, stable SPs cannot be found. 
Note that the (unstable) exact SP solution ® to the usual CCGL equation 
exists, unless both (32 and C2 vanish, not only at positive but also at negative 
values of the product /?ici ^D' , which is a drastic difference from the classical 
NLS solitons that exist solely in the case (3iCi > 0. Nevertheless, Fig. 4 shows 
that the parametrically driven CCGL supports, as well as its NLS counterpart, 
(stable) SPs only in the case (3ici > 0. 

2.3 Breathers 

As it was mentioned above, the SP solution suffers destabilization via the Hopf 
bifurcation at the border marked by rhombuses in Fig. 3, similar to what is 
known about solitons in the usual driven-NLS model ^ [3 |B| |Sj ■ The Hopf 
bifurcation gives rise to stable oscillatory pulses in the form of breathers. Similar 
breathing pulses are also found in the quintic CCGL equation JJj. While we did 
not aim to investigate the breathers in detail, we have found that they always 
keep a single-humped shape. With further variation of the parameters, the 
breather performing regular oscillations undergoes period-doubling bifurcations, 
which do not destabilize the pulse as a whose, but render its intrinsic oscillations 
chaotic. Figure 5 shows the time evolutions of the breather's amplitude, which is 
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(a) (b) (c) 
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Figure 5: Time evolutions of the breather's amphtude at (a) — C2 — 0.005, 
(b) /32 = C2 = 0.002 and (c) /32 = C2 = 0.001 in the case k = 4.5, ci = /3i = 
0.5,7= 1-5- 



equal to |u(L/2)|, for 7 = 1.5 and k = 4.5. The parameters (3i and ci are again 
fixed to be 0.5. The amplitude of the pulse exhibits a limit cycle oscillation 
at 132 — C2 — 0.005, a limit cycle oscillation with doubled period is seen at 
(32 = C2 — 0.002, and the time evolution becomes chaotic at (32 = C2 = 0.001. 



3 Periodic waves 

In the case when the condition Q does not hold and the zero solution is unsta- 
ble, SP solutions are irrelevant. However, it is natural to expect that the model 
supports spatially periodic stationary waves in this case, which may be quite 
relevant in the corresponding physical settings, provided that these waves are 
stable. 

Periodic solutions can be first analyzed by means of an analytical approx- 
imation, similar to what was done above for the SP solutions, i.e., adopting 
the ansatz ^ with constant 6. Then, an exact solution to the corresponding 
equation Q is 

u{x) = ^0 cn(ra, 9)e*^ , (17) 

where cn is the elliptic cosine with the modulus q ~ A^j \/ + and other 
constants are defined as follows: k = ^ (a^ + ^q)/2, = 2(fc -f 7Cos2(?) + 
^4(/c-F7COs26') -h2£;, and = -2{k + j cos29) + ^yA{k + j cos29) + 2E. In 
these expressions, E is a. free parameter of the family of the periodic solutions, 
which can be fixed if a certain period A is selected by periodic conditions. In 
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Figure 6: Comparison of the real and imaginary parts of the approximate an- 
alytical solution H17|l . i.e., AqCh (kx , q) cos 9 (solid curve) and AqCh (kx , q) sin 9 
(dashed curve), and a direct numerical solution, in the case with C2 = (32 — 0.02, 
k = —2,7 = 1.5, and the period A = 20/9. Actually, the analytical curves com- 
pletely overlap with the numerical ones. This solution is dynamically stable. 



terms of the above parameters, 

A = 4V2K{q)/^a^ + Al (18) 

where K{q) is the complete elliptic integral of the first kind. 

To determine the value of 9, we follow the same approach as in the case 
of SP, viz., we substitute the expression (|17|l into Eq. (|11|) . then multiply it 
again by the expression (|17|l . and integrate over the period (|18|l . The resultant 
equation was solved numerically, in order to find 9 for a given period A. 

In a direct numerical simulation, we have obtained a stable pattern with 
the period A = 20/9 at C2 = /32 = 0.02, fc = -2 and 7 = 1.5. In Fig. 6, we 
compare the direct numerical solution and the predicted approximate analytical 
solution for the periodic wave. A more general comparison of the analytical and 
numerical results is presented in Fig. 7, cf. a similar picture for the SP solutions 
in Fig. 2 [the average phase shown in Fig. 7(b) is defined by the expression 
pt)|) . in which the integration is limited to one period of the solution; in fact, 
the phase of the periodic solution is almost constant, unlike that of SPs]. The 
dependences in Fig. 7 terminate at points beyond which the numerical solutions 
could not be continued. As concerns the numerical solution, in the case shown 
in Fig. 7, a solution with A = 20/9 cannot be found for /32 > 0.04. Figure 
8 displays a time evolution at (32 = C2 — 0.045. The initial condition is the 
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Figure 7: Comparison of the amplitude (a) and average phase (b) of the sta- 
ble stationary periodic solutions, as found from the analytical approximation 
(continuous curves) and numerical solution (symbols). The dependences are 
obtained for the case f32 = C2, the other parameters being fixed, with k = —2, 
7 = 1.5, and the period A = 20/9. 



periodic solution with the period A = 20/9 for (32 = C2 = 0.02 shown in Fig. 6. 
The amplitude of the periodic solution with the period A = 20/9 decays to 
zero, and another periodic solution with period A = 20/6 appears as a stable 
solution. 



4 Conclusion 

In this work, we have introduced a model of the one-dimensional nonlinear 
dispersive medium with intrinsic loss, which are compensated by a parametric 
driving field. The model is a combination of the parametrically driven non- 
linear Schrodinger (NLS) and complex cubic Ginzburg-Landau equations, and 
has various physical applications (in particular, to optical systems). In the case 
when the zero background is stable, we developed an analytical approximation 
for a solitary-pulse (SP) solution, which was found to be in good agreement with 
direct numerical results. Unlike the well-studied case of the driven damped NLS 
equation, in the present case SPs have a nontrivial phase structure. Combining 
the analytical and numerical methods, we have identified the stability domain 
for the SP solutions. Generally, the increase of the diffusion and nonlinear-loss 
parameters leads to shrinkage of the stability region, but a nontrivial feature is 
that it may stretch in the direction of large values of the frequency mismatch. 
At one border of the stability domain, the SP is destabilized by the Hopf bifur- 
cation, which gives rise to a localized breather. Further period doublings make 
internal vibrations of the breather chaotic. In the case when the zero back- 
ground is unstable, hence SPs are irrelevant, we constructed stationary periodic 
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Figure 8: Time evolution of the real part of u in the case with C2 = /32 = 0.045, 
k = —2, 7 = 1.5. The period is A = 20/9 in the initial state, and the period is 
A = 20/6 in the final state. 



solutions, for which a very accurate analytical approximation was developed, 
and their stability was studied in direct simulations. 
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